Scaling near Quantum Chaos Border in Interacting Fermi Systems 
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The emergence of quantum chaos for interacting Fermi systems is investigated by numerical 
calculation of the level spacing distribution P(s) as a function of interaction strength U and the ex- 
citation energy e above the Fermi level. As U increases, P(s) undergoes a transition from Poissonian 
(nonchaotic) to Wigner-Dyson (chaotic) statistics and the transition is described by a single scaling 
parameter given by Z = {Ue a — uo)e 1 ^ 2 " , where uo is a constant. While the exponent a, which 
determines the global change of the chaos border, is indecisive within a broad range of 0.9 ~ 2.0, 
the finiteness of v, which comes from the increase of the Fock space size with e, suggests that the 
transition becomes sharp as e increases. 



PACS numbers: 05.45.-a, 05.30.Fk, 73.23.-b 
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Recently, the emergence of quantum chaos due to 
electron-electron interaction has attracted much atten- 
tion. While the subject has a longer history in nuclear 
physics, with the progress of modern nanofabrication 
techniques, it has entered condensed matter physics since 
a quantum dot system, for example, could be considered 
as an artificial atom with some physical parameters un- 
der control. Besides being of interest as its own right, the 
importance of the subject stems from the fact that it is 
related to a failure of the perturbative approach in inter- 
acting many-particle systems. Traditionally, the pertur- 
bative method has been one of the standard tools in theo- 
retical many-particle physics. The emergence of quantum 
chaos means a strong mixing of the unperturbed levels, 
thereby inducing breakdown of the perturbation series. 

A recent theoretical work by Altshuler et al. Q for the 
quasiparticle decay in a quantum dot has especially stim- 
ulated many theoretical investigations. In their paper [Q, 
the quasiparticle decay process was mapped to a single 
particle diffusion on the Bethe lattice, which is a non- 
perturbative treatment of the problem. They concluded 
that there is a transition to quantum chaos at a critical 
excitation energy e c ~ ^/T/A, where g is the dimensionless 
conductance and A is the mean level spacing between the 
single-particle levels. However, many authors afterwards 
have pointed out that the mapping to a Bethe lattice in 
ref. Q is too naive for a proper description and they ob- 
tained different results using other methods. The ongo- 
ing controversy could be summarised through two main 
questions, i) what is the relation between e and g at the 
quantum chaos border, and ii) whether the transition is 
sharp or not. Regarding the first question, there exist 
further different predictions such as e c ~ g 2 ^ 3 A (2jj|] and 
e c ~ 5 A/ In 5 As for the second question, Jacquod 
and Shepelyansky argued that the transition is smooth 
since the coupling between the Fock states is of nonlo- 
cal nature. The authors of ref. || concluded that the 
transition is smooth based on their numerical result for 
the local density of states and the participation number. 
In ref. ||, the same conclusion has been drawn by cal- 
culating the inverse participation ratio (IPR) for higher 



values of e by use of the, so called, layer model. On the 
other hand, Berkovits and Avishai M suggested a finite 
size scaling hypothesis, according to which the transition 
becomes sharp as e increases, which was based on their 
exact numerical results for small system size. Silvestrov 
jl| also proposed that the transition is sharp when the 
effective high-order interaction is taken into account. In 
addition, the authors of ref. J§| performed an analysis of 
their numerical result for the IPR and found that their 
results show features consistent with the prediction of 
ref. Q. However, their conclusion is not decisive enough 
concerning the sharpness of the transition. A reliable 
numerical test is, therefore, urgently needed to settle the 
issue. 

The main difficulty of numerical test is due to the fact 
that one should consider the regime of g 3> 1, which cor- 
responds to e > A. Since the size of the matrix to be 
diagonalized rapidly increases with e, one needs an al- 
ternative to the brute force method. The layer model, 
introduced by Georgeot and Shepelyansky , allows one 
to handle much larger system size (higher e) at a given 
computational cost by truncating out the Slater determi- 
nants which contribute little to a given eigenstate. 

In this paper, we calculate the level spacing distribu- 
tion P(s) for interacting fermionic systems up to e/A = 
27 by use of the layer model. The change of P(s) from 
the Poissonian to Wigner-Dyson statistics represents the 
transition from integrability to chaos. While our result 
does not allow us to resolve the controversy over the 
parametric relation for the quantum chaos border, i.e. 
question (i), it gives strong evidence that the transition 
becomes sharper as e increases. The finite size scaling 
(FSS) behaviour can be understood through a compari- 
son with an infinite dimensional Anderson model. 

Let us begin with rif spinless fermions on m single- 
particle levels for which the Hamiltonian TL is given by 
Ho + Tii with 



and Hi = £ U l j^ic\c s k c i Cj. (1) 

i<j,k<l 



H.q represents the noninteracting Hamiltonian. A single- 
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particle level energy with i = 1, 2, • • •, m is randomly 
chosen from the interval [(i — 1/2) A, (i + 1/2) A] with uni- 
form probability distribution. Uij.ki denotes the random 
two-body interaction matrix element and also has a box 
distribution [— U/2, 17/2]. When 7i is interpreted as the 
Hamiltonian for a quantum dot system, U is related to 
g as [/ ~ A/g 0JT^]. W is widely known as the Two- 
Body Random Interaction Model (TBRIM) and its full 
matrix size is given by N t = nf\/[m\(m — n/)!]. When TL 
is written in matrix form on the basis of the eigenstates 
of TLq, a diagonal element is the sum of the energies of 
occupied single-particle levels, while an off-diagonal cle- 
ment is nonzero only for two states which are different 
upto two particle-hole pairs. 

Since Nt increases very rapidly with n/ and m, di- 
rect numerical diagonalisation of TL is limited to rather 
small values of and m. Due to this difficulty, the 
layer model has been introduced by Georgeot and She- 
pelyansky (|] and the details are as follows. For each 
eigenstate of TLo, we define a sequence {fi} such that 
fi = when the level i is empty and fi — 1 when it 
is filled with a fermion. Then £/A = YliLi^fi ranges 
from n/(n/ + l)/2 to n/(m+ 1) — ?i/( n / + l)/2, which is 
approximately equal to the ground state energy and the 
highest excitation energy, respectively, in units of A. If 
we rewrite TL in ascending order of £ and set Uij,ki = 
for i + j 7^ k + I, the whole matrix reduces to a block 
diagonal form where £ is constant within each block. We 
call the submatrix with £/ A — n/(n/ + l)/2 = j the j-th 
layer model. In general, the matrix size of the j-th layer 
model, N(j), varies depending on n/ and m. However, 
for 1 <C rij <C m, N(j) for j < rif is determined solely 
by j and behaves as IO] 



N(j) ~ exp(ny/2j/3)/j. 



(2) 



The layer model hereafter is understood in such a sense. 
An eigenstate of TL with energy E can be written as a 
superposition of the eigenstates of Ho, the energies of 
which lie within the width T around E. T is much less 
than A, when U <C A, i.e. g> 1. This defines the valid 
regime of the layer model. 

One of the well established criteria for transition from 
integrability to chaos is the change of the level spacing 
distribution P(s) from the Poissonian P p (s) = e~ s to 
the Wigner surmise P w (s) = (7rs/2)e _7rs2 / 4 . To quantify 
the proximity of P(s) to either of the two, it is useful to 
define rj in terms of the variance of s as follows; 



n 



f o °° s 2(P( s )-P w ( s ))ds 
J °°sZ(P p (s)-P w (s))ds- 



(3) 



In this way, 77 is 1 for P p (s), for P w (s) and inbetween 
for an intermediate distribution. 

In our calculation, the Hamiltonian for the layer model 
for 19 < j < 27 has been constructed and numerical 
diagonalisation has been performed over 200-1000 disor- 
der configurations for each parameter set (J, g) . Corre- 
sponding matrix size ranges from 490 (j = 19) to 3010 



(j=27). To exclude the contribution of the tail states 
near the edge, 50 % of eigenvalues around the band cen- 
ter have been used to construct P{s). If we choose a 
smaller part of the eigenvalue set, e.g. 25 %, the result 
does not change significantly. 

The ?7's for 19 < j < 27 and for 1/200 < U < 1/10 
p"2[ are shown in Fig. 1. The gradual change from 
i] ~ 1 (Poissonian, integrable) to i] ~ (Wigner-Dyson, 
chaotic) with increase of U is found with all j. As j 
increases, the transition occurs at a smaller value of U, 
indicating that the interaction becomes more efficient for 
mixing the levels as the energy increases. The global de- 
pendence of the chaos border on the excitation energy 
may be found by plotting the data of Fig. 1 on a rescaled 
x-axis with x = Uj a [[L3|. If there exists a = ao such 
that all the data points lie in the same curve, the tran- 
sition is described by the relation e = XJ~ x l a ° and there 
is a smooth crossover instead of a sharp phase transition 
when e — * 00. However there does not exist such a value 
«o for our data. Instead, as shown in Fig. 2(a-c), there 
is a single crossing point where the curves meet one an- 
other and the slope at the crossing point becomes larger 
as j increases. This suggests that at sufficiently high en- 
ergy the transition to chaos is a sharp phase transition. 
In fact such a FSS feature with j is found for a broad 
range of 0.9 < a < 2. When a is tuned from 0.9 to 2, rj c 
(77 at the crossing point) decreases monotonically from 
0.8 to 0.2. For a < 0.9 or a > 2, the crossing point 
cannot be clearly identified. These results imply that 
while we cannot determine the exponent a which gov- 
erns the quantum chaos border p[ , the transition shows 
a FSS property regardless of the choice of a as long as 
0.2 <7 lc < 0.8. 

For a further analysis of this FSS behaviour, we as- 
sume that 77 is given by a function / of a single scaling 
variable Z such that 



r,(e,U)=f(Z) = f[(Ue a -u )e 



1/2 V] 



(4) 



where uq is a constant of the order of unity which is in- 
dependent of U and e. To understand the introduction of 
this scaling parameter Z, it would be helpful to compare 
the layer model with an Anderson Hamiltonian. One can 
think of the standard Anderson Hamiltonian defined on 
a graph, which has the same structure as the j-th layer 
model; there are N(j) vertices and a bond between two 
vertices exists if two corresponding Slater determinants 
of the j-th layer model have nonzero transition ampli- 
tude by the two-body interaction. The size of the graph 
can then be defined by the average minimum number of 
bonds between two vertices. It is given by ~ y/J, which 
is the same as the typical number of electron-hole pairs 
for the j-th layer model || . The spatial dimension of the 
graph in the limit of j — > 00 is given by 



\nN(j) 



(5) 
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We expect that there exists a transition to chaos as we 
decrease the disorder strength of the on-site energy for 
a fixed value of hopping energy between nearest neigh- 
bors. The transition will show a FSS behaviour with 
the system size L ~ j 1 / 2 and with v — 1/2, the corre- 
lation length exponent of infinite dimensional Anderson 
model [[TH). Returning back to the layer model, we expect 
that a similar FSS property exists since there also is the 
competition between the system size and a correlation 
length governing the chaotic property of the system. A 
difference between the layer models and the thus defined 
Anderson models is that the matrix elements are strongly 
correlated in the former, while they are completely ran- 
dom in the latter Jig] , Therefore, v of Eq. (4) is not 
necessarily equal to that of the infinite Anderson model, 
1/2. In passing, we note that Berkovits and Avishai R] 
introduced a scaling hypothesis very similar to Eq. (4) 
based on their numerical results for lower energy. 

According to Eq. (4), the slope of the curves r](U) for a 
given j behaves as ~ j a + 1 /(. 2v ) & ^ ^ ne crossing point and 
v is obtained by fitting this variation of the slopes for a 
given value of a. Fig. 3 shows the scaling plot obtained 
in this way with a = 1.5 (r) c = 0.56) and in this case 
v is given by 0.59. Such a procedure can be performed 
in the same way for 0.2 < r\ c < 0.8 (0.9 < a < 2) and 
the result for l/(2i/) is shown in Fig. 4. Taking account 
of uncertainties of all data, 1/(2^) lies between 0.5 and 
1.5 (0.3 < v < 1) over the whole range where the cross- 
ing point is identified. Therefore, though v also varies 
in a rather broad interval, our data evidently excludes 
the possibility of v being infinity, i.e. the possibility of a 
smooth crossover. One should note that this large varia- 
tion in v over a factor of 3 is the result of the fact that 
a is indecisive. If one can pin down the range of a in 
some other way than our method, or assume one of the 
previously proposed values (e.g. 1.5 ||,|| or 2 |l[]) to be 
valid, the uncertainty in v can be determined to a higher 
accuracy from Fig. 4. 

Finally, we discuss the sharpness of the transition when 
e is varied for a fixed value of U , since in real measure- 
ments with a quantum dot the conductance of the sample 
g ~ A/U is kept constant and the bias voltage, which cor- 
responds to e, is varied. We define the transition as sharp 
when the ratio of the transition interval Se to the transi- 
tion energy e c goes to zero as e c increases. From Eq. (4), 
e c is given by {u /U) 1 / a and Se is estimated using the 
slope of the curve r](e, U) as 



Se 



di] 
Ik 



(6) 



Therefore, (Se/e)\ t=ta ~ e^ 1 ^ 2 "' 1 and we find that the 
transition becomes sharp as long as v is not infinity. 
Therefore, the result of Fig. 4 shows that the transition 
to chaos takes place sharply if one uses a sufficiently clean 
sample, i.e. large <?, for a quantum dot. 

In summary, we have numerically investigated the 
emergence of quantum chaos in the interacting many- 



fermion systems using the layer model. While our result 
does not allow us to find the global dependence of the 
chaos border as a function of the interaction strength, 
we find that the transition becomes sharp at sufficiently 
high energy. This is also true when one observes the 
transition to chaos in a quantum dot as a function of 
bias voltage. 
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FIG. 1. rj as a function of the interaction strength U for 
the layer model with various values of j. 




(a) a =1.2 




FIG. 2. The same data as Fig. 1 with rescaled rc-axis of 
x = Uj" with a — 1.2, 1.5 and 1.8. A partial set of data from 
Fig. 1 is shown for clarity. 
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FIG. 3. Scaling plot of the data of Fig. 1 for a = 1.5. In 
this case 1/(2!/) is given by 0.85, i.e. v — 0.59. 




FIG. 4. 1/(2^) when r) c varies from 0.2 to 0.8. 



(c) a = 1.8 



10 12 14 



4 



